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We report calculations of the minimum energy pathways connecting 'CH 2 +N 2 to diazomethane and 
diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to 
HCN 2 +H, and of diazomethane to CH 2 N+N, The calculations use complete active space 
self-consistent field (CASSCF) derivative methods to characterize the stationary points and 
internally contracted configuration interaction (ICCI) to determine the energetics. The calculations 
suggest a potential new source of prompt NO from the reaction of ] CU 2 with N 2 to give diazirine, 
and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN 2 , which leads 
to HCN+N( 4 S) on the previously studied CH+N 2 surface. The calculations also predict accurate 0 
K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, 
respectively. © 1995 American Institute of Physics. 
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I. INTRODUCTION 

The reaction of singlet methylene ( ! CH 2 ) with N 2 along 
with the reaction of CH( 2 F1) with N 2 , which was previously 
studied by Walch, 1 are potentially important in the formation 
of “prompt” NO. 2 While triplet methylene is known to be 
important in some of the reactions leading to prompt NO 
formation, 2 the role of singlet methylene has not been con- 
sidered previously. There have been several previous studies 
of the CH 2 N 2 surface. An early study of the ground state and 
some of the low-lying excited states of diazomethane, at the 
ground state geometry, was carried out by Walch and 
Goddard. 3 More recent work has been carried out by Bold- 
yrev, Schleyer, Higgins, Thomson, and Kramarenko 4 
(BSHTK) who studied the diazirine and diazomethane 
minima and several saddle points on the CH 2 N 2 as well as 
the CHFN 2 and CF 2 N 2 potential energy surfaces. Guimon, 
Khayar, Gracian, Begtrup, and Pfister Guillouzo 5 (GKGBP), 
as part of a study on the decomposition of tetrazole, also 
reported a number of stationary points on the CH 2 N 2 surface, 
including a saddle point connecting diazirine and diaz- 
omethane. These studies were carried out at a lower level of 
theory than that used in the present study and did not provide 
sufficiently detailed information for combustion modelling 
studies. In this paper we report a detailed study of the path- 
ways for addition of ! CH 2 to N 2 to give diazomethane and 
diazirine and the pathway for rearrangement of diazirine to 
diazomethane. Calculations with extended basis sets are also 
reported which lead to accurate heats of formation for diaz- 
irine and diazomethane. 

In Sec. II we discuss qualitative features of the reactions 
considered here. Section ITT contains the technical details of 
the calculations, Sec. TV contains the results, and Sec. V 
concludes the paper. 

II. QUALITATIVE FEATURES 

Figure 1 shows the qualitative features of the orbitals for 
the reaction of 1 CH 2 with N 2 via a pathway in which the CH 2 
attacks the it bond of N 2 . In the ground state of meth- 
ylene ( ! CH 2 ) there is a substantial near degeneracy effect 


between the C 25 lone pair and the empty C 2 p orbital of a " 
symmetry. This leads to a pair of singlet coupled sp hybrid 
orbitals directed above and below the molecular plane, as 
indicated in Fig. 1(a). In the addition of ] CH 2 to N 2 , two 
bond pairs undergo major changes. These are the two elec- 
trons in the CH 2 C 2s pair discussed earlier and the in plane 
7r orbital of N 2 . By analogy to the reaction of CH( 2 II) with 
N 2 (Ref. 1) and the reaction of CH( 2 II) with H 2 (Ref. 6), the 
addition of 1 CH 2 to N 2 occurs via a pathway, as shown in 
Fig. 1(a), in which the CH and N 2 approach each other in a 



FIG. 1. The electronic structure of selected stationary points on the potential 
energy surface for CH 2 (U,)+N 2 : ( a ) jp, ; (b) min 1; (c) sp 3 ; (d) distorted 
diazomethane; (e) sp 2 ; (0 min 2. 
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TABLE I. Computed energies for stationary points on the surface. 


Geometry 

CIl a 

A£ b 





(a) c 



Reactants 

-148.427 2 1(- 0.469 40) 

0.0 



min I 

-148.471 56( 0.518 97) 

-24.0 



min 2 

-148.480 05<— 0.531 98) 

-33.9 



3 ch 2 +n 2 

-148.443 88(-0.484 81) 

-9.7 



Geometry 

CIi a 

A£ b 

CI2* 

A£ b 



(b) d 



Reactants 

-148.391 06(- 0.43026) 

0.0 

- 148.392 69(-0.431 14) 

0.0 

3 ch 2 +n 2 

-148.408 70(-0.44677) 

-10.4 



vdW 

- 148.390 96(-0.43381) 

-0.5 



sp 3 

- 148.398 69(-0.44458) 

-4.9 



sp I 

-148.383 37{ — 0.4291 8) 

4.7 

-148.385 66(-0.430 64) 

4.4 

min I 

-148.434 44(- 0.47866) 

-23.3 

-148.435 68(-0.479 61) 

-23.3 

spl 



- 148.359 65(-0.405 59) 

19.8 

min 2 

- 148.443 45(-0.49145) 

-32.9 


-33.2 

CH 2 N+N 



-148.291 37(-0.332 44) 

64.3 

hcn 2 +h 

— 148.282 39(— 0.32161) 

66.2 



Geometry 

CIl a 

A£ b 





(C) e 



Reactants 

-148.270 49(- 0.29892) 

0.0 



min ! 

-148.311 75(-0 34441) 

-21.4 



min 2 

— 148.319 34<— 0.35505) 

-29.7 



3 ch 2 +n 2 

— 148.290 72(— 0.31842) 

-12.3 



Geometry 

CASSCF 






(d) 



Reactants 

-147.998 52 

0.0 



sp 1 

-147.967 61 

23.5 



min 1 

-148.01541 

-3.5 



sp2 

-147.939 63 

40.8 



min 2 

-148.011 55 

-2.7 



vdW 

-148.004 30 

-1.9 



sp3 

-148.001 51 

2.2 




a The energies are in the form ICCI (ICCI+Q+148). 

Relative energies in kcai/mol including zero-point energy (see Table II). 


c cc-pVQZ basis set. 
d cc-pVTZ basis set. 
e cc-pVDZ basis set. 


nearly parallel arrangement. This non-least-motion pathway 
allows the formation of the diazirine molecule [Fig. 1(b)], 
while maintaining large overlaps between the components of 
the two bond pairs which change substantially during the 
course of the reaction (see Ref. 6). By contrast, addition of 
! CH 2 to N 2 via a C 2u pathway requires effectively breaking 
one bond and is a forbidden process. 7 

Another pathway leading to diazomethane involves ini- 
tial formation of a dative bond between one N 2s lone pair of 
N 2 and the empty C 2 p orbital of CH 2 [Fig. 1(c)]. This struc- 
ture evolves to diazomethane [Fig. 1(d)] by simultaneous 
donation of the N 2 s pair onto carbon and back donation of 
the C 2s pair into the tt system of diazomethane. This pro- 
cess involves no barrier. 

Figure 1 (e) shows the orbitals for the saddle point con- 
necting diazirine to diazomethane. The two radical orbitals in 
Fig. 1(e) arise from the CN cr bond of diazirine which is 
being broken. The doubly occupied orbital on the center N 


atom starts as a N 25-like lone pair, but becomes 5/?-like in 
the saddle-point region and eventually becomes a v orbital 
for diazomethane. The orbitals for diazomethane are shown 
in Fig. 1(f). From Fig. 1(f) it is seen that diazomethane has 
two 7 r electrons in a T orbitals and four rr electrons in a " 
orbitals (with respect to the molecule plane). The two singly 
occupied a " orbitals shown in Fig. 1(f) are singlet paired as 
a result of through bond coupling from the doubly occupied 
a ” orbital on the center N. 3 

III. COMPUTATIONAL DETAILS 

Several different basis sets were used in these calcula- 
tions. For the CASSCF derivative calculations, which were 
used to locate the stationary points, the polarized valence 
double zeta set of Dunning and Hay 8 was used. The basis 
sets for C and N are the (9s5p)/[3s2p) basis augmented by 
a single set of 3 d functions with exponent 0.75 for C and 


J. Chem. Phys., Vol. 103, No. 12, 22 September 1995 



Characterization of the minimum energy path for the reaction of singlet 
methylene with N 2 : The role of singlet methylene in prompt NO 

Stephen P. Walch 

Thermosciences Institute - MS 230-3 , NASA Ames Research Center, Moffett Field, California 94035 
(Received 11 October 1994; accepted 16 June 1995) 

We report calculations of the minimum energy pathways connecting , CH 2 +N 2 to diazomethane and 
diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to 
HCN 2 +H, and of diazomethane to CH 2 N+N. The calculations use complete active space 
self-consistent field (CASSCF) derivative methods to characterize the stationary points and 
internally contracted configuration interaction (ICCI) to determine the energetics. The calculations 
suggest a potential new source of prompt NO from the reaction of ! CH 2 with N 2 to give diazirine, 
and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN 2 , which leads 
to HCN+N( 4 S) on the previously studied CH+N 2 surface. The calculations also predict accurate 0 
K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, 
respectively. © 1995 American Institute of Physics. 


I. INTRODUCTION 

The reaction of singlet methylene ( ! CH 2 ) with N 2 along 
with the reaction of CH( 2 IT) with N 2 , which was previously 
studied by Walch, 1 are potentially important in the formation 
of “prompt” NO. 2 While triplet methylene is known to be 
important in some of the reactions leading to prompt NO 
formation, 2 the role of singlet methylene has not been con- 
sidered previously. There have been several previous studies 
of the CH 2 N 2 surface. An early study of the ground state and 
some of the low-lying excited states of diazomethane, at the 
ground state geometry, was carried out by Walch and 
Goddard. 3 More recent work has been carried out by Bold- 
yrev, Schleyer, Higgins, Thomson, and Kramarenko 4 
(BSHTK) who studied the diazirine and diazomethane 
minima and several saddle points on the CH 2 N 2 as well as 
the CHFN 2 and CF 2 N 2 potential energy surfaces. Guimon, 
Khayar, Gracian, Begtrup, and Pfister Guillouzo 5 (GKGBP), 
as part of a study on the decomposition of tetrazole, also 
reported a number of stationary points on the CH 2 N 2 surface, 
including a saddle point connecting diazirine and diaz- 
omethane. These studies were carried out at a lower level of 
theory than that used in the present study and did not provide 
sufficiently detailed information for combustion modelling 
studies. In this paper we report a detailed study of the path- 
ways for addition of ! CH 2 to N 2 to give diazomethane and 
diazirine and the pathway for rearrangement of diazirine to 
diazomethane. Calculations with extended basis sets are also 
reported which lead to accurate heats of formation for diaz- 
irine and diazomethane. 

In Sec. II we discuss qualitative features of the reactions 
considered here. Section HI contains the technical details of 
the calculations, Sec. TV contains the results, and Sec. V 
concludes the paper. 

II. QUALITATIVE FEATURES 

Figure 1 shows the qualitative features of the orbitals for 
the reaction of [ CH 2 with N 2 via a pathway in which the CH 2 
attacks the 7 r bond of N 2 . In the ] A } ground state of meth- 
ylene ( ! CH 2 ) there is a substantial near degeneracy effect 


between the C 2s lone pair and the empty C 2 p orbital of a" 
symmetry. This leads to a pair of singlet coupled sp hybrid 
orbitals directed above and below the molecular plane, as 
indicated in Fig. 1(a). In the addition of ] CH 2 to N 2 , two 
bond pairs undergo major changes. These are the two elec- 
trons in the CH 2 C 2s pair discussed earlier and the in plane 
7r orbital of N 2 . By analogy to the reaction of CH( 2 IT) with 
N 2 (Ref. 1) and the reaction of CH( 2 II) with H 2 (Ref. 6), the 
addition of ! CH 2 to N 2 occurs via a pathway, as shown in 
Fig. 1(a), in which the CH and N 2 approach each other in a 



FIG. 1. The electronic structure of selected stationary points on the potential 
energy surface for CH 2 (U 1 )+N 2 : (a) sp j ; (b) min 1; (c) sp 3 ; (d) distorted 
diazomethane; (e) sp 2 ; (9 min 2. 
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TABLE I. Computed energies for stationary points on the ! CH 2 +N 2 surface. 


Geometry 

CU a 

AE b 





(a) c 



Reactants 

-148.427 21 (-0.469 40) 

0.0 



min l 

-148.471 56(— 0.518 97) 

-24.0 



min 2 

- 148.480 05(— 0.531 98) 

-33.9 



j ch 2 +n 2 

-148.443 88(— 0.484 81) 

-9.7 



Geometry 

Cll a 

A£ b 

CI2 a 

AE b 



(b) d 



Reactants 

-148.391 06(- 0.43026) 

0.0 

- 148.392 69(-0.43 1 14) 

0.0 

’ch 2 +n 2 

-148.408 70(- 0.44677) 

-10.4 



vdW 

- 148.390 96(-0.43381) 

-0.5 



sp3 

-148.398 69(-0.44458) 

-4.9 



sp l 

-148.383 37(-0.42918) 

4.7 

-148.385 66(- 0.430 64) 

4.4 

min I 

-148.434 44(-0.47866) 

-23.3 

- 148.435 68(-0.479 61) 

-23.3 

sp2 



- 148.359 65(-0.405 59) 

19.8 

min 2 

-148.443 45(-0.49!45) 

-32.9 


-33.2 

CH 2 N+N 



-148.291 37(-0.332 44) 

64.3 

hcn 2 +h 

-148.282 39(— 0.32161) 

66.2 



Geometry 

CIl a 

A£ b 





(c) e 



Reactants 

-148.270 49(-0.29892) 

0.0 



min 1 

— 148.311 75(— 0.34441) 

-21.4 



min 2 

-148.319 34(-0.35505) 

-29.7 



3 ch 2 +n 2 

- 148.290 72(— 0.31842) 

-12.3 



Geometry 

CASSCF 






(d) 



Reactants 

- 147.998 52 

0.0 



sp 1 

-147.967 61 

23.5 



min 1 

-148.015 41 

-3.5 



sp2 

- 147.939 63 

40.8 



min 2 

-148.011 55 

-2.7 



vdW 

-148.004 30 

-1.9 



sp 3 

-148.001 51 

2.2 




a The energies are in the form ICCI (ICCI + Q+ 148). 

b Relative energies in kcal/mol including zero-point energy (see Table II). 


c cc-pVQZ basis set. 
d cc-pVTZ basis set. 
e cc-pVDZ basis set. 


nearly parallel arrangement. This non-least-motion pathway 
allows the formation of the diazirine molecule [Fig. 1(b)], 
while maintaining large overlaps between the components of 
the two bond pairs which change substantially during the 
course of the reaction (see Ref. 6). By contrast, addition of 
l CH 2 to N 2 via a C 2 „ pathway requires effectively breaking 
one bond and is a forbidden process. 7 

Another pathway leading to diazomethane involves ini- 
tial formation of a dative bond between one N 2s lone pair of 
N 2 and the empty C 2 p orbital of CH 2 [Fig. 1(c)]. This struc- 
ture evolves to diazomethane [Fig. 1(d)] by simultaneous 
donation of the N 2 s pair onto carbon and back donation of 
the C 25 pair into the it system of diazomethane. This pro- 
cess involves no barrier. 

Figure 1(e) shows the orbitals for the saddle point con- 
necting diazirine to di azomethane. The two radical orbitals in 
Fig. 1(e) arise from the CN a bond of diazirine which is 
being broken. The doubly occupied orbital on the center N 


atom starts as a N 25-like lone pair, but becomes sp - like in 
the saddle-point region and eventually becomes a it orbital 
for diazomethane. The orbitals for diazomethane are shown 
in Fig. 1(f). From Fig. 1(f) it is seen that diazomethane has 
two 7 r electrons in a ' orbitals and four 7r electrons in a” 
orbitals (with respect to the molecule plane). The two singly 
occupied a " orbitals shown in Fig. 1(f) are singlet paired as 
a result of through bond coupling from the doubly occupied 
a" orbital on the center N. 3 

III. COMPUTATIONAL DETAILS 

Several different basis sets were used in these calcula- 
tions. For the CASSCF derivative calculations, which were 
used to locate the stationary points, the polarized valence 
double zeta set of Dunning and Hay 8 was used. The basis 
sets for C and N are the (9s5p)l[ls2p] basis augmented by 
a single set of 3 d functions with exponent 0.75 for C and 
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TABLE n. Harmonic vibrational frequencies and zero-point energy for 'CH 2 +N 2 . Harmonic frequencies in 
cm obtained from CASSCF calculations with a pVDZ basis set. 



3214 

3100 

3307 

1951 

2329 

1660 

1461 

1509 

1546 

1081 

534 

1033 

477 

97 

1047 

72 li 

75 

836 

3348 

3176 

3436 

1009 

373 

1166 

389 

ZPE a 

80 

1018 


3313 

3360 

3253 

1595 

1902 

1659 

1491 

1495 

1438 

960 

1162 

1001 

747 

560 

3360 

888; 

328 

1021 

141 

3505 


3458 

1188 


1047 

419 



3432 

1570 

612 

928 

1251 

891 


0.029 45 0.025 68 0.034 28 0.029 05 0.03 1 7 1 

sp 3 


0.026 73 


0.019 76 


3142 

2342 

1471 

1180 

258 

248# 

3240 

1084 

232 


0.029 50 


’Zero-point energy in E H . The zero-point energy for the reactants is 0.022 97. 


0.80 for N (pVDZ). The H basis is ( 4 .y ) / [ 2 .s] augmented 
with a single set of 2 p functions with exponent 1.00. The 
basis sets used in ICCI calculations are the Dunning correla- 
tion consistent cc-pVDZ, cc-pVTZ, and cc-pVQZ basis sets. 9 

The calculations were carried out in C s symmetry with 
the mirror plane in the plane of the paper of Fig. 1. The 
active space for the CASSCF calculations included the eight 
electrons which are depicted in Fig. 1(a). This required six 
active orbitals of a' symmetry and two active orbitals of a" 
symmetry. The resulting CASSCF calculation had eight elec- 
trons distributed among eight orbitals. 

The CASSCF gradient calculations used the SIRIUS/ 
ABACUS system of programs, 10 while the ICCI calculations 
were carried out with molpro. 1112 All electrons were corre- 
lated except for the C 1 s and N 1 s like electrons with the 
restriction that no more than two electrons were allowed into 
certain weakly occupied orbitals. These orbitals were the 
highest three a and highest a" orbitals for the calculation 
denoted by CI1 and the highest two a' and highest a" orbit- 
als for the calculation denoted by CI2. A multireference ana- 
log of the Davidson’s correction 13 was added to the ICCI 
energies and is denoted by +Q. 

IV. DISCUSSION 

The computed ICCI energies obtained at the optimized 
CASSCF geometries are given in Tables 1(a), 1(b), and 1(c) 
for results obtained with the cc-pVQZ, cc-pVTZ, and cc- 
PVDZ basis sets, respectively. The zero-point effects were 
estimated as 1/2 the sum of the harmonic frequencies, which 
are given in Table II, and the zero-point effects are included 
in the relative energies, which are also given in Tables 1(a)- 


1(c). Thus, the relative energies are appropriate for compari- 
son to experimental results corrected to 0 K. The relative 
energies from Table I are also shown in Fig. 2. Table III 
gives the values of selected bond lengths and angles for the 
stationary points obtained here and also gives results from 



H2C— N=N min2 

FIG. 2. Computed energetics for the reaction of CH 2 ('A,) with N,. Note 
that while spZ is a saddle point at the CASSCF level, it is below the 
reactants energy at the ICCI level and it is probable that there is no barrier 
on this pathway. 
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HO. 3. Geometries for l 

per left sp , ; upper-nght min 1, center ten p 2 
left vdW; bottom-right sp y . 

omethane and 34.1 kcal/mol for diazirine. ^se values !^ 
elude zero-point effects and are appropriate for 0 K. Using 
heat of formation of 'CH 2 of 102.1 kcal/mol obtained from 
the heat of formation of 3 CH 2 of 93 kcal/mol and a smglet- 
tnplet splitting of 9.1 kcal/mol gives heats of formation of 
777 kcal/mol for diazirine and 68.0 kcal/mol for 
omethane. Laufer and Okobe 17 report the heat of formation 
oHiazinne as 60.6-66 kcal/mol. However there ,s an e~ 
her value of 79.3 kcal/mol by Paulett and Ettinger. The 
discrepancy between our calculated value and the value 
Laufer and Okobe is surprisingly large and our calculations 
support the value obtained by Paulett and Ettinger. au er 
and Okobe also report a heat of formation of 5 1 - ca mo 
for diazomethane. 19 This is also surprisingly far from our 
computed value. Calculations for the D e of N 2 using he 
same extrapolation method as used here give an extrapolated 
value of 226.8 kcal/mol compared to an expenmental va u 
of 228.4 kcal/mol. 20 As a further check on the computed 
results, we also computed the energy of 3 CH 2 +N 2 with each 
of the basis sets. The 3 CH 2 geometry was taken from the 
work of Bauschlicher, Langhoff, and Taylor who earned 

out a definitive study of the 

These results are also given in Tables 1(a) ( )• 


FIG. 5. Geometries 
omethane. 


along the reaction path connecting diazirine with diaz- 


seen that the cc-pVQZ basis set gives a ^parat.onof^ 
kcal/mol compared to 9.1 kcal/mol obtained in Ref. 21 with 
a comparable size ANO basis set. Taking these consider- 
ations together leads us to assign a maximum error of 2 
kcal/mol to our reported 0 K heats of formation for diazmne 

and diazomethane. . . . * 

From Ref. 20 it is expected that the barrier heights ob- 
tained with the cc-pVTZ basis set would change by less than 
0 5 kcal/mol if extrapolated to the basis set limit. This occurs 
because the barrier height involves an energy difference b 
tween a bond that is being broken and a bond that is forming 
and therefore there is a cancellation of errors^ Accordingly 
this extrapolation was not carried out for the barrier heig 
reported in Table 1(b). 

Table 1(d) gives energetics for the stationary points ob- 
tained at the CASSCF level with the pVDZ basis set As 
expected these energetics are much less reliable than those 
obtained from the ICCI calculations. The genetaltrends^re 
that the binding energies are underestimated and the ban** 
heights are overestimated. In spite of these obvious defers 
in the CASSCF energetics, comparison of CASSCF 
ICCI stationary points for the N+0 2 reaction s owe 
the CASSCF method does a very good job of predicting t e 
geometry in the directions perpendicular to the reaction . co- 
ordinate^ but that the saddle-point geometries can be shifted 
somewhat along the direction of the reaction coordinate. An 


Y 

H 


tttt 

FIG. 4. Geometries along the reaction path connecting 'CH 2+ N 2 with di- 
azirine. 







FIG. 6. Geometries along the reaction path connecting 'CH 2 +N 2 with d.a- 
zomethane. 
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extreme example of this effect, in the present calculations, 
occurs for sp 3, where the ICCI method finds an energy at 
the saddle-point geometry below that of the reactants. This 
has been interpreted to mean that the barrier is shifted far out 
into the entrance channel and in fact there is probably no 
barrier. In a case where the barrier persists at the ICCI level, 
the shift along the reaction coordinate is normally much less 
significant. 

We now return to consideration of other product chan- 
nels in the reaction of ! CH 2 with N 2 . Formation of CH 2 N+N 
is endothermic by 64 kcal/mol and involves a singlet-triplet 
crossing. Thus, this pathway is not a likely product channel. 
Dissociation of diazirine to CHN 2 +H can lead to either dou- 
blet or quartet CHN 2 . The quartet state ( q min 1 of Ref. 1) is 
17.0 kcal/mol above the doublet state (d min 2). However, 
CHN 2 + H is computed to be endothermic by 66 kcal/mol, 
with respect to reactants, for the doublet state of CHN 2 , and 
is also not a likely direct product. 

More likely processes leading to CHN 2 are hydrogen 
abstraction from diazirine or diazomethane, both of which 
are expected to be important products of the reaction of 
] CH 2 with N 2 . Hydrogen abstraction from diazomethane 
leads to an HCN 2 geometry like the dative structure of 
Manaa and Yarkony." 3 Previous theoretical studies 1 indicate 
that there is a large barrier between the dative structure and 
the region of the surface leading to HCN+N. Thus, this path- 
way is not expected to contribute to prompt NO formation. 
On the other hand, hydrogen abstraction from diazirine leads 
to an HCN 2 geometry like that involved in the doublet- 
quartet surface crossing region of the CH+N 2 surface. For 
example, the reaction 

diazirine + OH^CHN 2 + H 2 0 ( i ) 

is exothermic by 28.5 kcal/mol. In addition, assuming no 
vibrational relaxation, diazirine formed from 1 CH 2 +N 2 
would have about 38 kcal/mol of internal energy (energy of 
S P i minus energy of min 2), leading to a total of 66 kcal/mol 
of available energy. As the NN bond is considerably shorter 
in diazirine it is probable that CHN 2 will be formed with a 
large amount of energy in the NN stretching mode. [The 
quartet state of CHN 2 has a longer NN bond (2.377 A) than 
the doublet state (1.635 A), which in turn is longer than the 
NN bond in diazirine (1.253 A).] From Ref. 1 the saddle 
points leading to CH+N 2 on the doublet surface and to 
HCN+N on the quartet surface are 31 kcal/mol and 36 kcal / 
mol above d min 2, respectively. Thus, if this much energy is 
transfered into vibrational energy, these channels are ener- 
getically accessible. Furthermore, energy in the NN stretch- 
ing mode also favors the doublet to quartet crossing [In Eq. 
(1) the CHN2 is formed initially in the doublet state], leading 
to HCN+N( 4 S), as for the CH( 2 II)+N 2 reaction. 

V. CONCLUSIONS 

The reaction pathways for the reaction of ! CH 2 with N 2 
have been studied using complete active space self- 
consistent field (CASSCF) derivative methods to character- 
ize the stationary points and internally contacted configura- 
tion interaction (ICCI) to determine the energetics. 


! CH 2 +N 2 can lead to diazomethane, with no barrier, or to 
diazirine, with a barrier of 5 kcal/mol. Diazirine is computed 
to be about 24 kcal/mol below 1 CH 2 plus N 2 . Rearrangement 
to diazomethane involves a 20 kcal/mol barrier, with respect 
to reactants, and diazomethane is computed to be about 34 
kcal/mol below ! CH 2 plus N 2 . The CH 2 N+N and CHN 2 + H 
channels are computed to be endothermic by 64 kcal/mol and 
66 kcal/mol, respectively. Thus, these product channels are 
unlikely to be important. 

The computed potential energy surface would suggest 
that diazirine will be readily formed from *CH 2 plus N 2 . It is 
proposed that subsequent reaction of diazirine with hydrogen 
abstracters such as OH, which would be abundant in the 
combustion environment, would lead to doublet CHN 2 , with 
considerable vibrational excitation, especially in the NN 
stretch mode, which would facilitate doublet to quartet cross- 
ing, leading to HCN+N product. Thus, this pathway may be 
an important new source of “prompt” NO. 
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extreme example of this effect, in the present calculations, 
occurs for sp 3 , where the ICCI method finds an energy at 
the saddle-point geometry below that of the reactants. This 
has been interpreted to mean that the barrier is shifted far out 
into the entrance channel and in fact there is probably no 
barrier. In a case where the barrier persists at the ICCI level, 
the shift along the reaction coordinate is normally much less 
significant. 

We now return to consideration of other product chan- 
nels in the reaction of ! CH 2 with N 2 . Formation of CH 2 N+N 
is endothermic by 64 kcal/mol and involves a singlet-triplet 
crossing. Thus, this pathway is not a likely product channel. 
Dissociation of diazirine to CHN 2 +H can lead to either dou- 
blet or quartet CHN 2 . The quartet state (q min 1 of Ref. 1) is 
17.0 kcal/mol above the doublet state (d min 2). However, 
CHN 2 +H is computed to be endothermic by 66 kcal/mol, 
with respect to reactants, for the doublet state of CHN 2 , and 
is also not a likely direct product. 

More likely processes leading to CHN 2 are hydrogen 
abstraction from diazirine or diazomethane, both of which 
are expected to be important products of the reaction of 
*CH 2 with N 2 . Hydrogen abstraction from diazomethane 
leads to an HCN 2 geometry like the dative structure of 
Manaa and Yarkony. 23 Previous theoretical studies 1 indicate 
that there is a large barrier between the dative structure and 
the region of the surface leading to HCN+N. Thus, this path- 
way is not expected to contribute to prompt NO formation. 
On the other hand, hydrogen abstraction from diazirine leads 
to an HCN 2 geometry like that involved in the doublet- 
quartet surface crossing region of the CH+N 2 surface. For 
example, the reaction 

diazirine + OH— >CHN 2 + H 2 0 ( 1 ) 

is exothermic by 28.5 kcal/mol. In addition, assuming no 
vibrational relaxation, diazirine formed from l CH 2 +N 2 
would have about 38 kcal/mol of internal energy (energy of 
sp i minus energy of min 2), leading to a total of 66 kcal/mol 
of available energy. As the NN bond is considerably shorter 
in diazirine it is probable that CHN 2 will be formed with a 
large amount of energy in the NN stretching mode. [The 
quartet state of CHN 2 has a longer NN bond (2.377 A) than 
the doublet state (1.635 A), which in turn is longer than the 
NN bond in diazirine (1.253 A).] From Ref. 1 the saddle 
points leading to CH+N 2 on the doublet surface and to 
HCN+N on the quartet surface are 31 kcal/mol and 36 kcal/ 
mol above d min 2, respectively. Thus, if this much energy is 
transfered into vibrational energy, these channels are ener- 
getically accessible. Furthermore, energy in the NN stretch- 
ing mode also favors the doublet to quartet crossing [In Eq. 
(1) the CHN2 is formed initially in the doublet state], leading 
to HCN+N( 4 S), as for the CH( 2 n) + N 2 reaction. 

V. CONCLUSIONS 

The reaction pathways for the reaction of 1 CH 2 with N 2 
have been studied using complete active space self- 
consistent field (CASSCF) derivative methods to character- 
ize the stationary points and internally contacted configura- 
tion interaction (ICCI) to determine the energetics. 


1 CH 2 +N 2 can lead to diazomethane, with no barrier, or to 
diazirine, with a barrier of 5 kcal/mol. Diazirine is computed 
to be about 24 kcal/mol below ] CU 2 plus N 2 . Rearrangement 
to diazomethane involves a 20 kcal/mol barrier, with respect 
to reactants, and diazomethane is computed to be about 34 
kcal/mol below ! CH 2 plus N 2 . The CH 2 N+N and CHN 2 +H 
channels are computed to be endothermic by 64 kcal/mol and 
66 kcal/mol, respectively. Thus, these product channels are 
unlikely to be important. 

The computed potential energy surface would suggest 
that diazirine will be readily formed from ! CH 2 plus N 2 . It is 
proposed that subsequent reaction of diazirine with hydrogen 
abstracters such as OH, which would be abundant in the 
combustion environment, would lead to doublet CHN 2 , with 
considerable vibrational excitation, especially in the NN 
stretch mode, which would facilitate doublet to quartet cross- 
ing, leading to HCN+N product. Thus, this pathway may be 
an important new source of “prompt” NO. 
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